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I ' Abstract 

Q I We present a two-dimensional theoretical model for the slow chemical corrosion of 

O ' a thin film of a disordered solid by suitable etching solutions. This model explain 

different experimental results showing that the corrosion stops spontaneously in a 
k><( \ situation in which the concentration of the etchant is still finite while the corrosion 

;h ' surface develops clear fractal features. We show that these properties are strictly 

related to the percolation theory, and in particular to its behavior around the critical 
point. This task is accomplished both by a direct analysis in terms of a self-organized 
version of the Gradient Percolation model and by field theoretical arguments. 

Key words: Disordered solid, corrosion, gradient percolation, field theory, 
absorbing state phase transitions 



Preprint submitted to Elsevier Science 6 February 2008 



t=0 



t=l 



r,. 


r,2 


ri3 


'■,4 


riB 




r„ 


'■ia 


»-„ 


'-,4 


^is 


re 


^v 


ra 


rs 


r,o 


re 


'•y 


ra 


f-s 


r,o 


r^ 


^2 


^3 


r^ 


rs 


r^ < p(0) 
r^ < p(0) 
r, < P(0) 

Net(l)=Net(0)-3 


ri 




'•4 




V , Net(O) 

P(0)= ^^^^°^ 
V 




L 















Fig. 1. Sketch of the etching dynamics in a square lattice: the sites 2, 3, 5 are etched 
at the first time-step as their resistances are lower than p(0). At the same time the 
number of etchant particles in the solution decreases by 3 units, and a new part of 
the solid is uncovered. 

1 Introduction 



When an etching solution is put in contact with a disordered etchable solid, 
it corrodes the "weak" parts of the solid surface while the "hard" parts re- 
main un-corroded. During this process new regions of the solid (both hard 
and weak) are discovered coming into contact with the etching solution. If the 
volume of the solution is finite and the etchant is consumed in the chemical 
reaction, the etching power of the solution diminishes progressively and the 
corrosion rate slows down. When the solution is too weak to corrode any part 
of the hardened solid surface, the dynamics spontaneously stops. It is an ex- 
perimental observation (,l) that the etchant concentration at the arrest time 
is larger than zero. We show below that its value is strictly related to the 
percolation threshold of the considered solid lattice. We show also that the 
final connected solid surface has clear fractal features up to a certain charac- 
teristic scale a, i.e., the surface thickness. This is precisely the phenomenology 
that has been recently observed in experiments on pit corrosion of aluminum 
thin films (|l|). We show that the fractal features and the characteristic scale 
can be explained by the critical behavior of classical percolation around the 
percolation threshold (J2). 



2 The model 



A simple dynamical model, capturing the above mentioned phenomenology, 
has been recently proposed and studied (|3|;|j). The model is sketched in Fig. 1, 
and can be formulated as follows: 
(i) The solid is represented by a square lattice of sites with random resistances 



to corrosion rj uniformly distributed in the interval [0, 1]. It has a width L and 
a given fixed depth Y . The etching solution has a finite constant volume V , and 
contains an initial number N^ti^) of etchant molecules: the initial etchant con- 
centration is therefore C(0) = iVet(0)/V^. Experimentally the "etching power" 
p{t) of the solution at time t is roughly proportional to the concentration C{t). 
We assume p{t) = C{t) without loss of generality, and take pc < p(0) < 1, 
where pc is the percolation threshold of the lattice. 

(ii) At t = the solution is put in contact with the solid through the bot- 
tom boundary y = 0. At each time-step t, the solid sites belonging to the 
solid-solution interface with r^ < p{t) are removed, and a particle of etchant is 
consumed for each corroded site. Consequently, the concentration of the solu- 
tion decreases with t. Finally, depending on the connectivity criterion chosen 
for the lattice sites (e.g., first nearest neighbor connectivity for solid sites), 
m{t) new solid sites, previously in the bulk, come into contact with the solu- 
tion for the first time. At the next time-step, only these sites can be corroded, 
as the other surface sites have already resisted to etching and p{t) decreases 
with t. Calling n{t) the number of removed solid sites at time t, and N{t) the 
total number of removed sites up to time t, one can write 

/x nit) ,^^ Ei-on(t) ,„, N(t) 
pit + l)=p{t) - -^ = piO) - ^ii^-Ai =p(o) _ _^. (1) 

Note that, during the process, the solution can surround and detach finite solid 
islands from the bulk. The global solid surface is then composed by the union 
of the perimeter of the "infinite" solid, here called corrosion front, and the 
surfaces of the finite islands. At the end, the corrosion dynamics spontaneously 
stops at tf such that all the surface sites have resistances r > p{tf). 

The main features of the model, found through extensive numerical simula- 
tions (lattices with up to 2000 x 2000 solid sites in (j^), are: 
(i) The final value pj = p{tf) is slightly smaller than the percolation thresh- 
old Pc- The difference \pf — Pc\ ^ as {L/V)°''' (we take L/V ^ 1) with 
Op ~ 0.45, when the limits L,V ^ -|-oo are taken in the appropriate way (jj). 
All this can be explained through percolation theory which, in fact, implies 
that for p{t) < Pc the probability to have a connected path of solid sites all 
with r > p{t) (and then stopping corrosion) is equal to one in the infinite 
volume limit. 

(ii) The corrosion dynamics can be divided (see Fig. 2-a) into two regimes: 
a smooth regime when p{t) > Pc, and a critical regime when p{t) ~ Pc- The 
duration of the former is measured approximately by tc, defined by p{tc) = Pc] 
tc is found to scale with the ratio L/V in the following way tc ~ V/L. The 
duration of the latter, (tf —tc), scales as (V/L)"*/ with at^ — 0.55 (with a fur- 
ther linear dependence on log L of the scaling coefficient due to the extremal 
nature of t/ (|5|)). 
(iii) The surface of the solid in contact with the solution displays a peculiar 
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Fig. 2. a- Time decay of p{t) (continuous line), compared to the analytical approxi- 
mation p(t) = pq exp{—tL/V) in the smooth regime (dashed line) (4) in a triangular 
lattice. There is good agreement for p{t) > Pc- b- Three snapshots of the corroded 
solid (dark gray) by the solution (clear) at three different time-steps: initial, in- 
termediate and final. Black regions represent detached islands. The final corrosion 
front shows fractality up to a scale given by its thickness. 
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Fig. 3. a-Box counting measure of the fractal dimension of the corrosion front for 
length smaller than a. b-Final normalized histogram of the global solid surface: sim- 
ulations results (empty points) are compared with the analytical evaluation (con- 
tinuous line). 

roughening dynamics (see Fig. 2-b). In the first smooth regime, the corrosion 
has a clear mean direction given by the initial condition, and the corrosion 
front becomes progressively rougher and rougher, while finite detached islands 
are quite small. In the second critical regime, spatial correlations increase on 
the corrosion front and the dynamics generates a locally isotropic fractal ge- 
ometry, while the detached islands are larger. 

(iv) The final corrosion front is fractal with a fractal dimension Df = 1.753 ± 
0.05 (see Fig. 3-a) up to a characteristic scale a which is the final thickness of 



the front. On larger scales the corrosion front looks like a one- dimensional line. 
Also a scales with L/V as follows: a ~ [L/V)'"" with a„ ~ 0.57 ~ 1/-D/. All 
these properties are well explained in the framework of percolation theory as 
a self-organized version of the Gradient Percolation model (p;^). 
(v) Another important experimental observation captured by this simple model 
is the progressive hardening of the solid surface, due to the corrosion only of 
freshly discovered weak sites which, in turn, deplete the etchant concentra- 
tion in time permitting future corrosion of only weaker sites. This hardening 
effect is described quantitatively by the behavior of the normalized resistance 
histogram 0(r, t) of the global surface sites giving the density of such sites 
with resistance r. Obviously, 0(r, 0) = 1 with r G [0, 1], while (f){r,tf) is given 
in Fig. 3-b. As one can see, the final global surface is much more resistant to 
corrosion than the initial one. The time evolution and the final shape of 0(r, t) 
have been successfully obtained theoretically in 



3 Field theory approach and dynamical percolation 



In order to describe this model in the critical regime, it is possible to develop a 
phenomenological field theory approach "a la Landau" (see (|7|)) consisting in 
writing down the functional Langevin equation of the process around criticality 
directly from the analysis of the symmetries of the system. To this aim let us 
consider the following three different local densities or coarse-grained fields: 
(i) s(x, t) describing the local density in the point x at time t of material 
susceptible to be etched at any time after t (i.e., in the discrete model presented 
above, bulk solid sites and "fresh" solid surface sites freshly arrived to the 
solid-liquid interface). 

(ii) 5'(x, t) is the local density of passivated and inert material (i.e., surface 
solid sites having already resisted an etching trial and then immune or not- 
susceptible to be corroded at any future time-step), 
(iii) c(x, t) is the local density of corroded sites replaced by solution sites. 

The mean field equations (rate equations) describing the evolution of the av- 
eraged mean values of these magnitudes are to the leading order in the fields 



s{t) = -ac{t)s{t) 

g(t)=a(l-p(t))c(t)s(t) 

c{t)=ap{t)c{t)s{t) (2) 

where p{t) is the probability to etch an active site at time t, and a is a 
positive constant that we fix equal to one without loss of generality. The 
interpretation of the first equation is: in order for the density of susceptible 



sites to change (decrease) in a region, it is necessary to have locally both 
a non-vanishing density of etchant and raw solid material susceptible to be 
etched. This restricts the dynamics to active regions (i.e., part of the solution- 
solid interface) in which non- vanishing local densities of s and of c coexist. 
The second and the third relations in Eq. 2 say that an active site becomes 
either a c-site, with probability p{t), or alternatively, after healing, a g-site 
with complementary probability 1 —p{t). Note that, as s + c + g = 0, 
the total number of sites is conserved during the dynamics. We have written 
so far mean field equations in which spatial dependence and fluctuations are 
not taken under consideration. To this aim, it is convenient to introduce the 
activity field p(x, t) = c(x, t)s(x, t). From Eq. 2 it follows immediately that 

p(t) = -c(t)p(t)+p(t)s(t)p(t). (3) 



We now use for p{t) Eq. 1, noting that in terms of the coarse grained fields 
we have N{t) = /(ix[c(x, t) — c(x, 0)]. Using this observation, it is simple to 
derive the expression oi p{t) as a function of the activity field p(x, t): 



p{t) = p(0) exp 



1 * 1 r * 

— j dt' j d^p{-K, t') = p(0) exp - / dt'p{t' 

V J L 



(4) 



where p(t) = y J dn-pi^x., t) is the average activity at time t. Since p{t) is a 
function only of the average activity, it can be considered as a deterministic, 
smooth, positive and decreasing term in the final Langevin equation, even 
though p(x, t) is a stochastic field. Using Eqs. 2, 3 and 4 it is possible to write 
the mean field equation for the activity field: 

t 
dtp{^, t) = m{t)p{^, t) - p(x, t) J dt'{p{t) + p{t'))p{^, t') , (5) 





where m{t) = p(t)s(x, 0) — c(x, 0). In order to obtain the complete functional 
Langevin equation (i.e., the field theory) for the process in the neighborhood 
of the critical point (i.e., p(t) ~ pc), we have to add to Eq. 5 the noise term 
and the spatial coupling terms. The former is found simply by observing that, 
if in a local region there are k etchable sites in contact with the solution 
at time t, an average number p{t)k of them will be etched with a typical 



Poissonian fiuctuation of the order of Jp{t)M. This shows that the noise term 



is ~ J p{-K,t)ri{-K,t) (multiphcative noise), where ri{-K,t) is the typical white 
uncorrelated noise with zero mean and (r/(x, t)ri{x.', t')) = A6{'x — x.')6{t — t'). 
As a consequence of Eq. 5 and of the form of the noise term, it is possible to 
show that the only relevant term of spatial coupling in the neighborhood of 
the critical point is the diffusion term V^p(x, t). Since pit) is smooth, positive 



and, in the neighborhood of the critical point, p(t) ^ Pc > 0, at the end 
we obtain a field theory that belongs to the universality class of dynamical 
percolation (i.e., of a dynamical version of the classical percolation) (J9|): 

t 
dtpi^, t) = mp{^, t) -7p(x, t)j dt'pi^, t') + VV(x, t) + ^p(x, t)77(x, t) , (6) 



where 7 > 0, and m measures, in the dynamical percolation ,the distance 
p — Pc between the chosen constant occupation probability p and its critical 
value Pc- If r?7, > — active phase — the process generate an infinite cluster of 
occupied sites, and if m < — absorbing state (|lO( ) — the occupation process 
arrests exponentially fast (i.e., p(x, t) decreases rapidly to zero) leaving only 
a finite cluster of occupied sites. Since p(x, t) = is a solution of Eq. 6, once 
this state is reached the occupation dynamics (i.e., in our model the corrosion) 
stops spontaneously: p(x, t) = is a so-called an absorbing state, and at m = 
we have an "absorbing state phase transition". However, in our case (Eq. 5) m 
depend on time, and starting with a positive value, as p{t) decreases in time, 
it passes spontaneously from the active phase {m > 0) to the absorbing phase 
(m < 0) arresting rapidly after this crossover. The velocity of this crossover, 
that is, the duration of the critical regime, depends on the finite size of the 
system which in this way determine both the effective spatial gradient of p 
and the typical scale of fractality of the system, while the critical exponents 
are those of percolation theory (|2l). This shows, finally, that our model is a 
self-organized version of Gradient Percolation in any dimension. 
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